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ABSTRACT 

An analysis is presented of the assumptions and the applicability 
of the three theoretical methods for calculating condensations in ccsraic 
clouds where no pre-existing nuclei exist. The three procedures are; 
thermodynamic equilbrlum calculations, nucleation theory and a kinetic 
treatment which would take into account the characteristics of each 
individual collision. Thermodynamics provides detailed results on the 
composition temperature and composition of the condensate provided the 
system attains equilibrium. Because of the cosmic abundance mixture 
of elements, large super saturations in some cases and low pressures, 
equilibrium is not expected in astronomical clouds. Nucleation theory, 
a combination of thermodynamics and kinetics, has the limitations of 
each scheme. Kinetics, not requiring equilibrium, avoids nearly all 
the thermodynamics difficulties but requires detailed knowledge of 
many reactions which thermodynamics tvoids. It appears to be the only 
valid way to treat grain formation in space. To do this will require 
experiments on condensation that will provide an understanding of the 
process as it occurs in astronomical systems and also will provide the data 
on cross-sections, rates and related characteristics of the separate steps. 
A review of experimental studies concludes the paper. 
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INTRODUCTION 

Analyses of the cloud collapse phase of star formation, accumulation 
of solid grains into larger solar system objects, trapping of elements 
later incorporated into larger objects, thermal balance and temperature 
of clouds, are among the processes in which cosmic grains play an 
important role. The formation mechanism o. grains largely determine 
their composition, concentration as a function of time and position, 
radiative characteristics and accumulation probability. Consequently, 
a review of the status of formation mechanism of grains is very 
pertinent to this conference. 

In this paper I treat the case of condensation from a gas in uhe 
absence of any pre-existing grains or nuclei upon which condensation 
can occur. Some of the analysis is expected to apply to condensation 
on grains already present. The paper is not a summary of the 'state- 
of-the-art' in describing the formation of grains. It is, rather, an 
attempt to examine the basic assumptions of the various procedures used 
to describe condensation and to examine their applicability to astronomical 


systems. 



The case of lieturpenoouts nuclei! t ion on ions or of growth on pre- 
ejiistlitg grains is not discussed here. Iona may be important for siellfir 
or c li cumstellar condensation. Conclusions about the nature of the 
grain and the mechanism are ool: CKpected to be appreciably changed. 
Nucleatlon would occur at higher temperatures or lower concentrations 
than with a neutral vapor. This problem will be taken up at a later 
time. Very small grains may be present in some of the dense clouds 
considered here which formed from low density clouds. 

ir interstellar grains survive the collapse of a diffuse cloud, 
grains less than 100 S are probably present (Greenberg, 197J). Some 
aspects of the growth of gi'ains are given in Greenberg's report at the 
ronrerence. The process, in outline, is described by Greenberg (1973) and 
hv Field (1975) but detailed effects of gas-grain interaction (Hlrth 
mid Pound, 1963) are not included. As part of the Goddard study of 
cosmic grains, an analysis of this problem is intended. 

Tlu'oe roeehanlsma have been used to describe condensation. In 
historical order they are.* (1) thermodynamics, Wildt (1933), (2) kinetics, 
Lindhlud (1935), and (3) nudeation theory (Kamljo, 19i3f>, 1968; Donn et 
al. , 1968) which is n combination of thermodynamics and kinetics. These 
analyses have subsequently undergone considerable development. A com- 
prehensive review of the thermodynamic procedure has been given by 
Grossman and Larimer (1974) . A brief review of thermochemistry and 
nudeation theory has recentlv been published by Salpeter (1977). 

Kinetic analyses rc-qulre the most detailed input data and few 
general treatments (ter Hoar, 1944; KramcM- and ter llaar, 1946; 
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van de Hulst, 19A9; Hoyle and Wickramaslngle, 1962; Bonn, 1976) have 
been attempted. The papers listed indicate the complexity and con- 
sequent limitations of a kinetic analysis. 

IIIERMODYNAMIC EQUILIBRATION CALCULATIONS 

The simplicity of thermodynamics, in which the equilibrium composi- 
tion and phases in a system depend only on the final state, allows 
detailed results to be calculated. A determination of the chemical 
composition of the gas and any condensates require knowledge only of 
the free energies of all species which may occur. Procedures for 
obtaining the data and carrying through the calculations are given in 

the review of Grossman and Larimer (1974). They adopted a cloud of 

“4 

"solar" composition at a pressure of 10 atmospheres. Figure 1, taken 
from that paper, summarizes their condensation computations. The chart 
shows Che fraction of each element that has condensed and the mineral 
in which is occurs, as a function of temperature. Of the three methods 
for treating condensation only thermodynamics permits such detailed and 
comprehensive results to be obtained. Indeed, for most elements, it is 
the only procedure that yields any results at the present time. 

Unfortunately, there is a price that must be paid for this informa- 
tion; Che system must be in thermodynamic equilibrium. To the degree 
that the system is out of equilibrium, calculated results will deviate 
from the actual values. Several factors exist which may lead to dis- 
equilibrium. The most important among these appear to be the low 
pressures in clouds, and formation of metastable species which may 
persist Indefinitely. The effect of each of these factors will be 
)iiTc*nttinl'i'cl by the non-ex Istonrp In Che gns of the nominal moleoulp 
III I ti<« roiiitoimeil mimoiil. 
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Condensation Is the tendency of the system to attain et|uillbrluin. 

As with many chemical reactions an activation bsTier eKlsta which may 
cause a non-equilibrium state to be "frozen-in". With vapor phase 
transitions, the requirement for supersaturation makes condensation 
an inherently non- 'jqullibr ium process which can have aiynlf-cant effects 
(Blander and Katz, 1968; Blander and Abdul Gawad, 1969, ^l.tnder, 1971). 
The experimentally determined supersaturation for iron (Frurip ,*fnd Bauer, 
1977) Is 20 at 2000 K and Increases to 500 at 1650 K, The theoretically 
and experimentally well established need for supersa duration thus 
Introduces a major deviation from equilibrium. For condensation of iron 
at 1600 K, the iron vapor has Co be undercooled by 300*^. The kinetic 
effects of low pressure and possible metastable species will be con- 
sidered in Section IV where an analysis of the kinetics of condensation 
is treated. It is to be emphasized that classical thermodynamics 
predicts the equilibrium composition and phase distribution and not 
the results of a phase transition which, for the reasons Indicated, 
may yield non-equilibrium, metastable products. 

III. NUCLEATION THEORY 

llie essence of nucleation theory is the transition from the gaseous 
to the 'ondensed phase, liquid or solid. In its usual formulation 
(Frenkel, 1955; Hirth and Pound, 1962; Abraham, 1974) it is a combina- 
tion of thermodynamics and kinetics. The process can be represented 
by the set of equations: 

+A^ t J = 1, 2 jA (1) 
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ORIGINAL I'Aui .. 

OP POOR QUAliltV 

where in the monomer, A^, a cluster of j molecules or atoms and Aj* 
the critical cluster or condensation nucleus . The distribution n (Aj ) , 
for j not too near j*, is given by thermodynamics. Monomeric addition 
of an existing gaseous molecule is generally assumed, although, in some 
cases, e.g,, carbon, disuers and trimera may be involved. Only one type 
of condensible niolcule occurs and a cluster changes size only by 
colllsional capture or evaporative loss. Several treatments of binary 
systems have recently been made (Kiang et al., 1973; Shugard et al,, 
197A ; Wilemski, 1975). Kiang and Cadle (1975) treated ternary 'ystems. 
In these analyses multicomponent particles form without reaction. 
Condensation with reactions to form grains for which the molecular 
composition is different from that of the vapor has been considered 
explicitly by Courtney (1964) and in a more general but formal sense 
by Hlrschfelder (1974). Characteristics of clusters which occur in 
these two investigations are not knowi\ preventing application to 
actual systems. 

Salpeter (1974) proposed that polymers incident on a growing 
cluster may fragment through release of the binding energy of the 
particle. Part is captured and the escaping fragment carries away 


the excess energy, stabilizing the now larger cl-.’Ster. This may be 


an Important suggestion as we shall see later. 

For the fonnation of the most stable major minerals e.g., AI 2 OJ, 
CaTji^O^ (Grossman and Larimer, 1974) , nucleatlon is expected to be 
seriously inhibited. The building blocks blocks are small molecules 
e.g., AI 2 O, TiO, CaO, not the larger, nominal molecule of the crystal. 
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The nominal molecules yielding very stable macrosccpic crystals c^i,.iot 
exist until clusters of many molecules occur to form the bulk material 
possibly In a very poor crystalline array. 

At the low pressures in astronaralcal systems two additional factors 

-4 

enter. At a total pressure of 10 atmos, the partial pressure of 

-9 

condensible species will be below 10 atmos. The Interval between 
auccesslv collisions of condensible species with a small cluster will 

be greater than one second whereas the Interval for all molecular 

-5 

collisions will be about 10 s, Tlierefore, the rate of cluster growth 
is slow on molecular time scales and the stability of small clusters 
becomes an important factor. The binding energy per atom for metal 
clusters has been studied theoretically by several investigators and 
one experimental measurement for iron exists (Freund and Bauer, 1977). 

The results are summarized in Figure 2 where it is seen that the binding 
energies of very small clusters is about 20% of that for the bulk 
metal. No data exists for oxides. The stability and hence concentra- 
tion of these clusters will be much less than would be found using the 
bulk binding energies. 

A sec< ad important effecu of low pressure Is a marked deviation 
from tnermal equilibrium. The consequences for chemical reactions 
In the interstellar medium was pointed out uy Donn (1973). An 
an.ilysis of non-equilibrium molecular vibrational energy in late 
type giant and supergiant stars was discussed by Thompson (1973) . His 
cri-Lerion for a non-equllibriuni distribution, radiative relaxation rate 
equal to colllsional rate, is qualitative and very conservative. However, 
.JO win use the .same simple measure for deviation from equilibrium here. 
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At p » 10 acroos, the molecular colliaion frequency Is 10 b . The 

efficiency of translational to vibrational energy exchange is very 

poor (Cottrell and McGoubrey, 1961; Kondratev, 1964a), for sample 

molecules at 1000*^K an effic' ncy of 10 ^ may be adopted, i.e,, 10^ 

collisions are required for energy transfer to occur. This yields a 

transition probability for collisional excitation of about 100 s 

Representative radiative vibrational transition probabilities are in 

the range 100-1000/sec. mol. (Penner, 1959). Therefore, significant 

deviations from a Boltzmann equilibrium distribution in the sense of a 

-4 

deficiency in excited states will appear at 10 atmospheres. 

In summary, nucleation theory makes the following assumptions: 

(1) pressures are high enough for thermodynamic equilibrium to exist and for the 

-5 

time scale of nucleation to be short, /v 10 s, (2) there are unrcoctlve, 
condensible species, (3) the condensed molecule exists in the vapor, 

(4) no reactive species exist. As we have showi, none of these assump- 
tions apply to astronomical systems. 

Even with the simplified approximations, nucleation theory has only 
been extensively used for graphite (Kamljo, 1966; Donn et al., 1968; 

Walker, 1975, Tabak et al . , 1975). Mineral condensation sequences based 
on undercooling required by the large supersaturations have been derived 
by Blander and his associates in the references given earlier. 

IV. ‘dT NETIC MECHANISMS 

The well known distinction between thermodynamics and kinetic 
theory shows up in the astronomical condensation problem. Whereas 
thermodynamics can make detailed calculations without needing detailed 
processes and data as long as equilbrlum exists, kin* tics requires no 
special conditions to carry through its calculations. The price to do this 


Is a complete knowledge of oil Import, nC processes, energies, cross-- 
sections tind other relevant data. The consequence Is that calculations 
for specific condensates have been utteiapted only for ice luixtures 
(van de Hulst, 1949) and graphite (Hoyle and Wickraiusingle, 1962). 

An anaiyais which examined for the first tljne the formation kinetics 
of diatomic molecules as the initial step in grain formation was 
given by ter Hoar (1942) and Kramers and ter Hanr (1946) . An amorphous 
particle consisting of essentially all elements was proposed. Cernuschi 
(1947) reconsidered some of the basic steps in grain formation after 
the initial clusters formed. In 1962, Hoyle and Mickramasingle discussed 
graphite formation in the atmosphere of cool carbon stars. They made 
use of exjjerimental values of the heat of reaction of small carbon 
polymers and theoretical value for the equilibrium vapor pressure of 
carbon. 

Wliat we wish to emphasize here, with regard to these investigations, 
is the difficulty, initially, of correctly writing down the series of 
proce.sses that occur. Following this, the serious problem of obtaining 
energies, cross-sections and rates required for the calculation must 
be resolved. In the papers mentioned, mechanisms were highly simplified. 
In spite of the large concentration of reactive molecules, only those 
steps were included which added to the cluster, e.g., for graphite 
forjiiatlon, only reactions of the form: 

C + C t c (2) 

n n+1 

were included. A much more frequent step will be 

k 

H + C -+• Products 
n 
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(3) 


wltere chis products, CH, CH^ etc. (Bulooch and Olnndet, 197b) are 
produced at a rate k. Anorher expected reaction oeries ia: 


CH + C *>• C + H 
m n n+1 R> 


( 4 ) 


where may be an unatable fragment which yieldn H + H^. Thlti is the 
type ot process bv which pyrolytic grapltlte if produced. The poBsible 
significant effects of liydrocarbon reactions on carbon condensation was 
mentioned by Donn et ai . , . They neglected these mechanisms Sie- 

cause of the lack of knowledge of how to treat them. Similar considera- 
tions apply to the analysis of ter Hear (19<^i2) and the other kinetic studies 
Several other factors having important consequences on the kinetics 
of grain formation are briefly noted here. At low pressures three-body 
collisions which are needed to stabiliae the association of small 

moiecules are very Infrequent. The ratio of two-to-three body collisions 

“22 

is approximately 10 n (Kondratev, 196f»b) where n is the total 

3 -4 “8 

particle density per cm . In the pressure interval 10 to JO atreosphtres 

“7 -11 

this ratio drops from 10 to 10 . Salpeter s suggestion of fragmenta- 

tion of the incident molecule can be effective, the escaping fragment 
carries away the excess binding energy and allows part of the incident 
molecule to add to a eJuster^with less than about eight atoms. 

A second factor that may enhance condensation is the low vibrationaJ 
temperature at low pressures which depletes the h gher vibrational 
states of the molecules. Tlic primary energy mode affecting condensa- 
tion appears to be the Internal energy of the molecule, cluster or 
surface rather than the molecular kinetic energy (McCarroll and Ehrlich, 
1964; Shade, 1*564). Tlius, at high kinetic temperatures, the nen-Bol t .u’iiann 


dlstribucion raujine that tha aystani ia affectively much cold<-*r. Thla 
may be the duniluant factor allowing gralna to form in "hoe" clouds. 

A almllar effect applicable to amall aggregates which behave uk bulk 
solids was pointed out by Arrhenius and De (1973). 

Finally, a major contributor to condensation may be Ionization 
which Is wall known to be very effective in causing nucleatlon In many 
systems (Hrrth and Pound , 1963) . A more general possible consequence 
of ionization on condensation It* clouds has been considered by Arrhenius 
and Alfven (1970; rue also LeJmert (1970)). Those papers consider 
ginln formation where plasma effects dominate clouds processes whereas 
this review is concerned with regions when ionization has limited, 
specific consequences. 

V. REMARKS ON THE iHEQRETICAL MECHANISMS 

The analysis of the three mechanism of investigating condensation 
is now summarized. Thermodynamics yields a precise, detailed evolution 
of condensation in clouds if the system were always at equilibrium. 
However, because of the fact that condensation Is a transition process, 
not an equilibrium state, various factor may Interfere with the 
equilibrium state being attained. Of particular significance are the 
low pressures characteristic of astronomical systems, and the large 
supersaturation required to drive the condensation process. Nucleatlon 
theory also suffers from the requirement for thermodynamic equilibrium. 
In addition, it has been developed only for very simple laboratory 
systems. The basic assumptions are not valid for astronomical cases. 

A generalization of nucleation theory, more valid for present circum- 
stances may be made (Bonn, 1976). This procedure will require much 





more detailed knowledge than can readily be obtained and leade ub to 
the third ucheiiie - a kinetic treatment. Kinetics Is valid under any 
circumstances if all the steps and quantities are known. For n low 
pressure, cosmic abundance dywtem, we are very far from beiiu; able to 
carry out a detailed calcnlation. A simple, extreme calculation in which 
all refractorv atoms and jaolecules add to the cluster yields a boundary 
condition grain. This was done previoualv (Donn, 197b) wleh the resultant 
coiBpoBltion Fe^Mgj^raAlNl (f'lO) ^ . Tliia procedure and result Is somewhat 
analogous to that of ter Haar (1974.!) who was dealitig wit!i low Cetsperai ure 
diffu : ’uds where the total cone en-.. rat ion was about unity. The 
Ce.Upo8itlon for the grain lust given may be compared with the thermo- 
dynamic ’“esults of Figure 1. Because of the considerable deviations 
from thermodynamic equll Ibr luw, those c.ilculations will not be accurate 
descriptions of the grains. 

The general character of the thermodynamic results are ex- 
pected to provide guidelines loi the grain coiytpositon and condensa- 
tion temperatures. Inliililt ton of nucle.ition will muse calculated 
temperatures to be upi-et liieits. l-'I ament s , e.g., aluminum, titanium 
and calcium, are not expec'‘t‘d lo h.ive distinct condensation tempera- 
tures or condense In unique i rystalllne minerals. The tende.icv may 
be for groups of elc-j«.nC8 to condense together over a ruange of temper. i- 
tures in an apior pluaia aggrcg.ite. Solid .state diffusion at high temiperu- 
tures could cau>-'e cinneal I'ing, le.iding to a niul t verystnlline grain (Stephen 
1978). The actual gr.iin composl t ion can be expected tp fall between 
the two extremes. As pointed out in the earlier paper (Donn, 197b> a 
unique cc'iiposltion ie unlikely. ft Is to be emphasized that the nature 
of the system and the factors that favor oondensatlon indicate nen- 


equlllbrium solids will form. 



VI. CONDENSATION EXPERIMENTS 


It should be clear that condensation experiments appropriate to the 
conditions in astronomical systems are needed. Experiments on condense- 
tion of smokes have b'ien carried out by Lefevre (1967 , 1970) who produced 
grains of iron, iron oxide, silicon carbide and silica by on arc discharge 
in argon. Kamijo et al., (1973) used heaters, arcs, and plasma jets to 
varporlae a variety of similar materials, again in argon, b’teyer et al. 
(197^) used an arc to vaporize and recondense silica. In these experiments , 
condensation takes place under uncontrolled conditions with unknown 
temperatures and pressures. Except for iron oxide and silicon carbide, 
there is only a single condensible species. The analysis of condensation 
conditions in clouds in the previous sections has shown that these 
e jeriments are valuable for producing small grains of various materials 
and studying their properties but they do not contribute much to our 
understanding of Tstronomical condensation. Meyer (1971) has experi- 
mented with condensation on surfaces initiated by sputtering of selected 
targets by argon atoms. His procedure and results are difficult to 
relate to smoke formation in clouds. Of significance to the present 
analysis is Meyer's conclusion that some films formed by sputtering 
magnesium silicate targets could not be identified with any known 
magnesium silicate compound. Perhaps the experiment most pertinent 
to grain formation In astronomical clouds was Che attempt by Dot f eld 
and Hudson (1973) to nucleate graphite. They simulated a cool stellar 
atmosphere with an expanding beam of a hydrogen-methane mixture 
equilibrated at temperatures on the order of 2000 K but obtained a negative 
result v,*th only an upper limit of 10 percent of carbon condensing. 
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Turther work along this line could yield data on the kinetics and 

products of carbon condensation In liydrogen rich atmospheres. 

Currently, there are two series of experiments on nucleating smokes 

from vuporlaed refractory materials which extend the work already referred 

to in the direction of being moiC realistic for cosmic clouds. Thu worlc 

of Arnold and his colleagues were reported at the symposium of Stephens 

and is described elsewliere in these proceedings. To date their experiments 

have been done at pressures above an atmosphere. At Goddard, Kenrick 

Day and I vaporized magnesium and silicon monoxide in separate, adjacent 

crucibles (Day and Donn, 1978). These were independently heated, allowing 

different Mg/Si ratios to be obtained In the vapor. An ambient argon 

atmosphere at about 2 Torr pressure could be heated to control the con- 

o 

densation temperature. Temperatures from 0 to 500 C were employed. 

The strong dependence of condensate yield on temperature, described 

below, indicates that the ambient temperature was the controlling 

parameter. Calculated partial pressure of the reactant Mg and SiO 

vapors were between 0.1 and 1 Torr. The smoke yield decreased ns the 

ambient temperature was raised. At 500°G no smoke was collected unless 

the crucible temperatures were increased substantially, raising the 

reactant partial pressures. At 500*^C the supersaturation for experiments 

5 

in which no grains were collected is estimated to be greater than 10 . 

At all temperatures the grain size was a few hundred angstroms. The 
particles were amorphous, with variable Mg/Sl ratios. Their infrared 
spectra sho. ed structureless 10 and 20 |ira features very similar to those 
found in astronomical sources. When the grain., were annealed in vacuum 
at lOOO^C for 1 hour, a well developed forsterite (Mg^SiO^) spectrum 
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developed. Figure 3 ahowe the InL’rared apeetrum oE comlemjcd and 
annealed gruina. The dependence of smoke yield on temperature, Che 
high auperea turn cion and Che amorphous, non-acoiclviometric character 
of the grains, all support the kinetic mechanism for grain fornintion. 
Those expeviujents are being extended to Include ocher metals and more 
complex mixtures as well as ionization and other phenomena. The experi- 
ments will also be made more quantitative. 

The concept of organic grains pi'oposed by Moyle and Wickramasinghe 
(1977) has not been considered here. The kinetic treatment is expected 
to apply to any composition. 

Tlie ambitious, ultl*,fiate goa. of our investigation of grain fonmation 
at Goddard la to develop a method to pei-mit approximate but realistic 
calculations of conde" atlon in astronomical systems. 
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FIGURE CAPTIONS 


Figure 1. Summary of thermodynamic equilibrium calculations from a 

nebula of solar composition, P = 10 ^ atinos. (from Grossman 
and Larimer, 197^0 • 

Figure 2. Binding energy of clusters normalized to bulk solids 
(from Anderson, J. Chem. Phys. jM, 4046, 1976). 

Figure 3. Infrared spectrum of smoke condensed from SiO + Mg mixture 

at 300 K. initial condensote; heated to 800 K 

in vacuum; .... heated to 1300 K in vacuum. 
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